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ABSTRACT 

We study the formation of molecular hydrogen, after the epoch of re-ionization, 
in the context of canonical galaxy formation theory due to hierarchical clustering. 
There is an initial epoch of H2 production in the gas phase through the H~ route 
which ends at a redshift of order unity. We assume that the fundamental units in 
the gas phase of protogalaxies during this epoch are similar to diffuse clouds found 
in our own galaxy and we restrict our attention to protogalactic disks although 
some of our analysis applies to multi-phase halo gas. Giant molecular clouds are not 
formed until lower redshifts. Star formation in the protogalactic disks can become 
self-regulated. The process responsible for the feedback is heating of the gas by the 
internal stellar radiation field which can dominate the background radiation field at 
various epochs. If the gas is heated to above 2000 — 3000-ftT the hydrogen molecules are 
collisionally dissociated and we assume that in their absence the star formation process 
is strongly suppressed due to insufficient cooling. As we demonstrate by analysis of 
phase diagrams, the H2-induced cool phase disappears. A priori, the cool phase with 
molecular hydrogen cooling can only achieve temperatures > 300K. Consequently, 
it is possible to define a maximum star formation rate during this epoch. Plausible 
estimates give a rate of 0.2 — 2M & yr~ l for condensations corresponding to la and 
2a initial density fluctuations. For more massive structures, this limit is relaxed and 
in agreement with observations of high redshift galaxies. Therefore, the production 
of metals and dust proceeds slowly in this phase. This moderate epoch is terminated 
by a phase transition to a cold dense and warm neutral/ionized medium once the 
metals and dust have increased to a level Z « 0.03 — 0.1Z Q . Then: (1) atoms and 
molecules such as C, O and CO become abundant and cool the gas to below 300K ; 
(2) the dust abundance has become sufficiently high to allow shielding of the molecular 
gas and; (3) molecular hydrogen formation can occur rapidly on grain surfaces. This 
phase transition occurs at a redshift of approximately 1.5, with a fiducial range of 
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1.2 < z < 2, and initiates the rapid formation of molecular species, giant molecular 
clouds and stars. The delayed initiation of the cold phase in the interstellar medium 
of protostellar disks at a metallicity of Z <; 0.1Z Q is consequently a plausible physical 
reason that the formation phase of the stellar disks of the bulk of the galaxies only 
occurs at a redshift of order unity. The combination of feedback and a phase transition 
provides a natural resolution of the G-dwarf problem. 

subject headings: cosmology: theory - galaxies: formation - galaxies: evolution - 
molecular processes 



-4 - 



1. Introduction 

Despite numerous searches for molecules at high redshift, the positive results are sparse. The 
only confirmed observations of molecules at high redshift seen in emission are associated with 
two gravitational lenses: FSC 10214 + 4724 (Scoville et al. 1995, Frayer 1995) and the Cloverleaf 
(Barvainis et al. 1995). In absorption, H2 has been seen in PKS 0528 - 250 (Foltz et al. 1988) at 
a column density of 10 18 cm~ 2 , but in no other damped Lyman Alpha system. This is surprising 
since the column densities of atomic hydrogen inferred for damped Lyman Alpha systems are 
~ 10 21 cm -2 . In comparison with diffuse clouds in our own galaxy such columns seem sufficient 
to produce H2. Galactic diffuse clouds cover the disk with a covering factor between 20-50% 
where the lower limit is from the Copernicus data (Savage et al. 1977) and the upper limit comes 
from assuming all the infrared Cirrus is associated with diffuse clouds. Other molecular species 
including HCO+, HCN, CO, C 18 0, HNC, CN, CS and H 2 CO have been seen below redshift ~ 1 
in the millimeter band (Combes & Wiklind 1996). These absorption techniques are effective for 
any redshift, provided a bright millimeter continuum background exists. Nevertheless, all the 
absorbers known to date are in lensed systems with correspondingly small impact parameters. 

At lower redshift z < 0.2 there are well established observations of many molecular species in 
emission in the nuclei of starburst galaxies (c.f Scoville & Soifer 1991). For large nearby galaxies 
such as M31 and NGC 891 (c.f. Young 1990) CO maps have been obtained. For our own galaxy, 
about a hundred different molecular species have been detected in interstellar space. Many di- 
and tri-atomic molecules can be observed toward diffuse clouds through the visual and millimeter 
absorption lines they produce in the continua of nearby stars. Some quasars (~ 30%) used as 
millimeter callibration sources are known to exhibit absorption lines as well (de Geus et al. 1996). 

These data pose interesting questions about the prevailing chemical balance and the physical 
structure of the ambient medium as a function of redshift (c.f. Norman & Braun 1996). Previous 
work on this subject has focused on very high redshift in the epoch prior to reionization where 
light-element hydrides such as DH and LiH have been discussed in the context of observations of 
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deuterium abundances and faint structure in the cosmic background radiation field (Maoli et al. 
1996, Dalgarno et al. 1996, Stancil et al. 1996). The importance of H 2 cooling for the formation of 
the first stars has been discussed by many authors (c.f. Haiman et al. 1996, Tegmark et al. 1996). 
Although the self-consistent formation of H 2 above a redshift of z ~ 6 due to gas phase processes 
yields very low abundances of molecular hydrogen H 2 /H ~ 10~ 6 (Lepp & Shull 1984), these may 
be sufficient to sustain the collapse of initial perturbations. In the post re-ionization epoch, the 
problem of H 2 formation has been discussed by Black et al. (1987) for the QSO PHL 957 — a 
damped Lyman Alpha system with an atomic HI column density of 2.5 x 10 21 cm ~ 2 — where the 
H 2 abundance was more than five orders of magnitude lower than values representative of our own 
galaxy (c.f. Levshakov Sz Varshalovich 1985, Lanzetta et al. 1989). Their analysis indicated that 
two essential parameters are: (1) the background radiation field which is inferred to be an order of 
magnitude higher than typical intergalactic values at z = 2.3, and (2) the dust-to- gas ratio which 
is constrained to be no more than ~ 1/3 of typical galactic values. Many new results have recently 
become available on the UV background radiation field (cf. Haardt & Madau 1995), the dust 
content of damped Lyman Alpha systems (Pettini et al. 1994, Fall & Pei 1994) and the merging 
history of dark matter halos (Kauffmann h White 1993). QSO absorption line studies have also 
made remarkable progress in the last few years. The high redshift (z > 2) dependence of the HI 
column density distribution function contains valuable information on the early stages of structure 
formation during which large systems are likely to be in the process of aggregating themselves 
from smaller sub-units. The obvious candidates to focus on are those systems with HI column 
densities above N(HI) > 10 21 cm~ 2 (Wolfe 1995, Storrie-Lombardi et al. 1994). These systems 
are often associated with protostellar disks and are believed to contain a reasonable fraction of the 
baryons of the universe. 

These new developments motivate us to re-investigate the formation of molecular hydrogen 
after re-ionization and to assess its role in the standard galaxy formation theory due to hierarchical 
clustering and the subsequent formation of stars. If the latter process is to occur, a diffuse 
cool (T <; 500.fr) phase needs to be supported which ultimately leads to the formation of dense 
molecular clouds and stars. Obviously, the formation of stars produces an internal radiation 
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field which will heat its surroundings and influence the chemical balance. For star formation 
to continue, the cool phase should persist. In fact, this process will determine the maximally 
sustained star formation rate. If the associated time scale is long (comparable to a Hubble time), 
the abundance of metals and dust will increase slowly until a threshold is reached at which the 
cooling is dominated by atoms and ions like O and C + , as well as molecular species like CO. In 
addition, H2 formation on grains proceeds more rapidly and dominates in regions where the bulk 
of the material is in neutral form. Also, the increased columns of dust absorb the radiation of 
newly formed stars and sustain a star formation cycle as in our own galaxy. The onset of the 
bulk of the star formation at low redshifts and moderate metallicity is consistent with current 
constraints on disk formation. In this paper we have restricted our attention to protogalactic disks 
since they are the likely sites of the cold molecular phase we are focusing on. Multi-phase halos 
have been discussed extensively for galaxy formation (c.f. Kang et al 1991, Ikeuchi & Norman 
1991) and our results in sections 2.5, 3, 4, and 5 are relevant for such studies although we do not 
pursue them further here. 

This paper is organized as follows. In Section 2, the formation and subsequent evolution of 
structure in a CDM-dominated Universe is discussed, and the redshift dependence of the UV 
background radiation field, the stellar radiation field, the gas pressure, the radiation pressure and 
the ionization parameter are presented. The gas-phase epoch of molecular hydrogen formation is 
calculated in Section 3 including the production of metals and dust. The feedback effect due to 
star formation is analyzed in Section 4 and upper limits to the star formation rate and metal and 
dust production rate are derived. In section 5 we discuss the transition to a multi-phase ISM and 
the onset of rapid star formation. The epoch when H2 formation on grains becomes important 
is discussed in Section 6. The formation of metals and subsequent molecular cooling to massive 
molecular clouds is also discussed there. Section 7 discusses the implications for the study of the 
interstellar medium in protogalaxies and star formation, as well as the G-dwarf problem. We 
associate the onset of star formation in disks with both a redshift, z ~ 1 — 2, and a metallicity, 
Z ~ 0.03 — 0.1 as also indicated in a perceptive paper by Wyse and Gilmore (1988). Where 
possible we give reasonable analytical estimates and then use the full (numerical) details of the 
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recently calculated metagalactic and stellar radiation fields in our calculations. 

2. Basic Cosmological Model 

2.1. Galaxy Formation 

For the present investigation the specific details of structure formation are of minor importance 
and we adopt the semi-analytical relations presented by White & Frenk (1991) which agree well 
with more detailed N-body simulations (cf. Lacy & Cole 1993, 1994; Kauffmann et al. 1993, Heyl 
et al. 1995). In the numerical calculations, we also include the merging history of dark matter 
halos through the conditional probability that a halo of mass Mo at redshift zq has previously 
been in a halo of mass M\ at z\, as presented in Kauffmann & White (1993). 

Using simple top-hat models, the collapse to virial equilibrium of a perturbation that has 
become non-linear at redshift z results in a density enhancement 5. In CDM scenarios the fiducial 
value is S = 178 (Narayan & White 1988). We assume that for protogalactic disks to form, cooling 
takes place and a fraction of the gas originally contained in the mass perturbation collapses to 
a centrifugally supported disk with a collapse factor of A -1 , where A is the spin parameter from 
tidal torques (cf. Peebles 1993) with a canonical value of A = 0.07. Note that there are details of 
the cooling and infall to the protogalactic disk that we have omitted but they are not essential to 
the argument given here. For a density perturbation of mass M, the resulting column density is 
given by 

»-(f)' t ' ,, ^»* + "'' (1) 

yielding 

2 1 

^ixW^V ( \ V(l + ^) 2 cm- 2 , (2) 

vo.oiy V5 x io n M y v ' ' v ' 

with Qb t g the baryonic mass fraction in the protogalactic disk relative to the total perturbation. 
Our numerical estimates are sensitive to this parameter. In fact, £l^ g depends on the merging 
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history of the dark matter halos and is a function of redshift. Kauffmann (1996, her Figure 12) 
presents the redshift distribution of the baryonic fraction including the conditional probability for 
halo merging quoted above. There is a clear maximum between redshift 2 and 3 for a biasing 
parameter 1 < b < 2. A value of b = 1 — 1.5 seems most consistent with the latest derivations 
of Storrie-Lombardi & MacMahon (1996). Therefore, we adopt her results in the numerical work 
presented below. 

To characterize the density, we assume that the scale height, H, to disk size, R, denoted by 
rf = H/R, of the collapsed systems is approximately constant with a typical value of rf ~ 1/100. It 
follows that the density in the disk is approximately equal to 

n = (JA-V^bW 1 + z f ( 3 ) 

yielding 

- 5 (v)(^) (1+3)3 cm "' (4) 

The canonical mass value is chosen to be M = 5 x 10 11 M Q . Our estimates of the redshift 
of molecule formation are only weakly dependent on this choice if we focus on the bulk of the 
matter which is contained in the la and 2a density perturbations, M < 10 12 M Q . We also assume 
that at high redshift the formation of stars has not been effective in locking up the baryons. The 
ratio tp = (gas/gas + stars) should be of order unity for these proto-galactic objects. Of course, 
in more evolved systems such as our Galaxy the ratio is only ~ 10%. The actual frequency 
distribution of objects of mass M is determined by the cosmological model in the context of the 
Press-Schechter hierarchical clustering picture (c.f. White & Frenk 1991). For a cosmic density 
parameter Q = 1, it is assumed here that the total baryonic fraction in galaxies fi& )ff , at a given 
redshift, cools to a disk while the dark matter remains in the halo. The fraction of baryons in 
galaxies is typically 10% of the total number of baryons available. The cosmic density of baryons 
ribQ is found to be n^o = (3Hq /8irG)(Qb/ /imp) = 10 _6 (Ob/0.1)/i 2 , where Hq is the Hubble constant, 
h = Ho/100kms~ 1 Mpc~ 1 , Clj, is the baryonic fraction of matter in the Universe and \x denotes the 
reduced mass of the primordial gas for a helium abundance by number of 10%. 
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The above formulation is for single objects of a given mass and it is therefore worthwhile to 
formulate a statistical criterium which is more directly related to the initial conditions of structure 
formation. Equation (1) can also be written as 

N = n r 5 2 / 3 \- 2 n big (l + z) 2 , (5) 

where ro is the initial radius and r$ = (3M/47T) 1 / 3 . Using the Press-Schechter formalism (see 
White & Frenk 1991, Equations (1) and (2)), the abundance of halos n(M, z) with mass between 
M and M + dM is given by 

n(M ' z)= y 2( M^ )y dioiM exp( -T ) (6) 

where y = <5 C (1 + z)/a(M), we take 5 C = 1.68 and for standard CDM cosmology 
ct(M) = 16.36~ 1 (1 - OSdO^ 1 + 0.4814rg- 2 )- 10 , where r is in units of Mpc and is given 
by r = l(M/lO 12 M )(fi/i 2 )~ 2 / 3 Mpc. The function a(M) can be fitted by an approximate power 
law a = o"o(M/Mo) a , leading to power law distribution functions n(M, z) oc M~ a ~ 2 (l + z) for the 
halo mass distribution below a characteristic turnover mass defined by a(M*(z)) = S c (l + z). Gas 
cooling and resultant disk formation can only occur if the cooling time of the virialized gas in 
the dark halo at a temperature T = 2 x 1O 6 (M/1O 12 M )(1 + z)K is shorter than the dynamical 
time at a given redshift. This leads to an effective cooling mass as a function of redshift denoted 
M coo i(z). We show in Figure 1 the halo abundance distribution n(M, z). Also shown in the lower 
M — z plane of Figure 1 are the turnover mass M* and the limiting mass of the dark halo, M coo i, 
in which there is gas, with Qb,h = 0.05, that can cool at the given redshift (see also Section 5). 



2.2. The Background UV Radiation Field 



The background UV radiation field is well approximated for 0.5 < z < 2.5 by 

J= J- 21 (l + zf+i 



(7) 
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where J_2i is the fiducial value for the intergalactic background in units of 1CT 21 erg cm~ 2 s _1 
sr _1 Hz -1 . The value of q lies between 0.5 < q < 1, and reflects the uncertainty in the detailed 
spectral shape of the UV continuum of quasars. Current observations seem to favor a value 
q = 0.5 (Tytler et al. 1995). For z < 3, Equation (5) is accurate to within 30% over the 912-2000 
A range which is the relevant wavelength region for the H2 photodissociation process. Above a 
redshift of three, the quasar population turns over and the metagalactic background decreases 
with increasing redshift. 

In order to more clearly distinguish the effects of H2 dissociation and ionization heating, the 
background radiation field is separated into two components, one above and one below the Lyman 
limit at 912^4, denoted by J > and J < , respectively. The actual calculation of the radiation field 
as a function of both wavelength and redshift requires a numerical integration of the cosmological 
radiative transfer equation in the absorption line forest (Haardt & Madau 1996). We call this 
radiation field J5 which is the one used in the numerical computations. 



2.3. The Stellar Radiation Field 



When star formation is initiated in the protogalaxy, the internally generated radiation field 
can become significant. For example, at 1000 A the Draine-Habing estimate for the radiation field 
in our Galaxy, Jd-h, gives an effective radiation field of J-2i,D-H = 30 — 50. We calculate the 
luminosity, L*, and spectral energy distribution for a continuous star formation rate of lM yr _1 , 
with a Salpeter Initial Mass Function (IMF), and a lower mass cut-off of IM Q (Leitherer & 
Heckman 1995). To obtain an intensity, we assume that the stars are being formed uniformly 
troughout a disk with area A whose nominal value is 100/epc 2 . In the numerical results we include 
the more detailed estimates of Kauffmann (1996, her Figures 8 and 13) for the redshift-dependent 
size of the disk and the radial HI column density distribution. We denote the stellar radiation 
field by J* which is given by 

J* = U (?) , (8) 
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where S is the actual star formation rate. The total radiation field is then 



J = Jb + </*• 



(9) 



Figure 2 presents the shape and magnitude of the overall radiation field J for various epochs 



pronounced for wavelengths longward of the Lyman limit, leading to variations in the relative 
contributions to J < and «/>. 



The visual extinction to distant QSOs for damped Lyman Alpha systems is not more than 
a few tenths of a magnitude (Fall h, Pei 1994). Fundamental limits have been obtained on the 
depletion of elements like Zn and Cr in damped Lyman Alpha systems at a redshift of order ~ 3. 
Inferred dust-to-gas ratios, £ g d, are of the order of ^ g d ~ 0.01 — 0.1£o where £o is the mean galactic 
value (Pettini et al. 1995). The metallicities at these high redshifts are down by approximately 
the same factor, Z ~ 1/30Z Q . An accurate analytical model for the redshift dependence of ^ g d 
and Z is difficult to construct due to the essential role played by star formation. Nevertheless, 
this dependence should be included since it is a crucial ingredient in our analysis of molecule 
formation with cosmic epoch. In our numerical work, the metallicity is calculated from the star 
formation rate in the closed box limit and also in the limit of the ratio of stars to gas being small 
i.e. l — if> = star/(star + gas) ~ star/gas << 1. Consequently we can write 



where y « 0.02 is the yield of metals like C and O (Woosley & Weaver 1995). 

The production of dust follows by assuming that a fraction, £j, of metals is instantly locked up 
into grains (typically silicate cores). Therefore, we find that the dust-to-gas ratio is given by 



and various star formation rates. For S > 0.02 the stellar component dominates which strongly 
increases the intensity in the 912-1110A region where H2 is dissociated. The increase is most 



2.4. Dust and Metallicity Content 




(10) 




(11) 
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As we shall show later, the critical metallicities in our analysis turn out to be ~ 0.03 — 0.1Z Q . 
Therefore, we have neglected the significantly lower initial metallicity that may have come from 
the formation of the halo. 



2.5. Gas Pressure, Radiation Pressure and Ionization Parameter 

The pressure, P, in the mid-plane of the proto-galactic disk can be estimated by assuming 
that the mass in the disk is in the gas phase and using the Archimedean formula P = pgh where g 
is the gravitational acceleration normal to the disk. For a thin disk g is given by g = 27rGS where 
G is Newton's constant and £ is the surface density of matter in the disk. Therefore, one can write 

P = 2irn 2 m 2 p GN 2 , (12) 

which can be written using equation (1) as 

P = 2-h-%7r l iG^mj5^\- 4 n b M%nl{l + z) 4 , (13) 

with numerical value 

P=Y2m(^A ( \ V(l + ^) 4 cm~ 3 K. (14) 

vo.oiy V5xio n M y v ' y 1 

where P = P/k, and k is Boltzmann's constant. The (1 + z) A dependence of the pressure is slightly 
misleading due to the essential role played by £lb,g- I n f ac t, due to the maximum in the redshift 
dependence of the latter, the sizes of protogalactic disks are smaller at redshifts z > 3 compared 
to z ~ 1. Note also that at redshift z ~ 1, the pressure calculated from the weight of the gas is an 
order of magnitude greater than the thermal pressure calculated from taking nT. In fact, if one 
takes the radial dependence of the column density distribution into account and averages over only 
the inner 3 kpc of the disk, then the gravitational pressure is two orders of magnitude larger than 
the thermal pressure. This is also the case in our Galaxy and consequently there must be some 
turbulent pressure that holds up the gas (Boulares & Cox 1990, Norman h Ferrara 1996). The 
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turbulent pressure is due to kinetic motions from supernovae and superbubbles ultimately due to 
massive star formation. For low star formation rates the turbulent pressure may not be effective 
and the mid-plane pressure acting on the thermal component may increase leading to a density 
increase of more than an order of magnitude. This needs to be kept in mind when discussing the 
phase diagrams describing the state of the protogalactic gas. 

The pressure, P^, in a virialized halo with a gas fraction is given by 

A = 80f ) f fg^V(l + .) 2 cm-K. (15) 



,5xlO n M y V0.05. 

We therefore neglect the halo pressure relative to the disk pressure. More detailed models require 
the consideration of the pressures and density profiles in the halo and the role of cooling flows to 
precipitate material from halo to disk. These complex issues will not be addressed further here 
(c.f. Norman and Meiksin 1996) 



The radiation pressure is given by 



Prad = (16) 



and we verify that it is always less than the gas pressure in the proto-disks for each of our 
calculations. We define the ionization parameter used here by U = J/n with n the total number 
density, and a numerical value defined by U-21 = J-21/n. The dependence of the pressures and 
the ionization parameter are shown in Figure 3. The radiation pressure is typically much less than 
the gas pressure for the epochs studied here, but is of the same order of magnitude for z ~ 0. The 
ionization parameter varies by a factor of five over the 1 < z < 3 range, because of the turn-over 
in the quasar population around z = 2.5. 



3. Gas-Phase Formation of H2 with Background Radiation 

We now turn to examine the epoch when molecular hydrogen is produced in the gas phase. 
We first use a simple model (Donahue Sz Shull 1991) to find the approximate analytical expressions 
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for the molecular hydrogen abundance through the gas phase H route. 

'" UA = 1.4 x 10- 5 T a88 [/_2ix, (17) 



n 



where x = n e /n is the electron abundance and is given by 

x = 0.02T 0A2 Ul 21 . (18) 

These equations assume that the H~ abundance is in equilibrium and its dominant destruction 
channel is the reaction with neutral hydrogen yielding H2. Conversely, the dominant destruction 
channels of molecular hydrogen are UV photodissociation and collisional dissociation at 
temperatures T > 2500 K. A simple model for the ionization heating yields 

T = 3 x W 3 U°-lfn - 12 , (19) 

and the molecular hydrogen abundance is given by 

= lO-^ffn " 14 . (20) 

For our two component model, where the radiation field is split into a component longward of 
912A and one shortward of the Lyman limit, we find the temperature to be 

T = 3 x 10 3 (C/< 21 ) a27 ([/> 21 ) a54 n ai2 ^, (21) 

and the molecular hydrogen abundance to be equal to 

' TL ' ] 10- 2 (U< 21 f 85 (U> 21 )- - 3 n 015 . (22) 



n(H) 

In the numerical work, the optically thin H 2 dissociation rate assumed in the above equation has 
been corrected for the fact that at abundances of the order of lcolumn densities in excess of 10 19 
cm~ 2 , H 2 self-shielding plays a role (see below). 



4. Gas-Phase Formation of H 2 with Star Formation 
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4.1. Feedback 



It is generally expected that feedback mechanisms may be important in galaxy formation (c.f. 
Norman & Ikeuchi 1996). The following model is considered here. If the internal star formation 
generates a dominant radiation field then the gas will be heated to above 2000 — 3000i<C and the 
molecules will be collisionally dissociated. We assume that the absence of molecular hydrogen 
formation would strongly limit the efficiency of star formation. From the simple model used above 
we find that for the stellar radiation field the temperature can be estimated to be 

giving a limiting star formation rate of 

^'^To^' 125 ^ 11 *!^ M ° yr "'' <24> 

Using the estimate of the density from equation (2) we find that 



- °- 5 <io^> 125 (v) L1 {ml (looy (1 + z) " M ° yr ~' (25) 



For a Schmidt star formation law where S oc N p with, 1 < p < 2, this corresponds to a critical 
column density 



l.i /n. \ i.i N 



^-^K^ 1 -^) (M T (1+ ^ cm " 2 (26> 

For a protogalaxy of a given mass we find that, using the arguments outlined in Section 2, the 
radius of the disk decreases with redshift as 

«=ior^^W^y J (i+,)- kP c (27> 



5xlO n M y V0.01 



yielding, with A = ttR 2 , 



• 1000.fr \v) 5xlO n M J V0.01 



For star formation to proceed at all, S should be a factor of a few smaller. 



4.2. Bursts of Star Formation 



We have analyzed the feedback effects in terms of a fixed, continuous star formation rate. 
We now discuss the burst mode where we envisage the following cycle: the star formation rate is 
sufficiently high, so that all molecules are dissociated; the star formation turns off; the molecules 
reform; and another burst occurs. 

The time scale, td, for the dissociation of molecular hydrogen, assuming no self-shielding, is 

t d = 2 x 10 5 JZ21 yr (29) 
and the timescale, tf, for the formation of molecular hydrogen is 

^ = 4x10 *(— ) {—) yr - (30> 

Consequently, the duty cycle, D = (tf/td), of on-burst to off-burst is 

D = 8 x 10~ 3 C/- 2 i( T ) (31) 

Since the amplitude of the radiation field is set by the level of star formation rate it follows that 
in any burst the mass of the stars formed M*^ urst is constant and given by 

M* Mrst = St d w 1O 6 M . (32) 

Thus, small bursts can occur but with a very low mean effective star formation rate. It should 
be emphasized that once a multi-phase medium has been established, bursts can be sustained 
(Spaans & Norman 1996). 



5. The Multi-Phase Structure of the ISM 



5.1. Numerical Model 
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We adopt the heating and cooling curves presented by Donahue & Shull (1991). We include 
cooling by [01] 63 fim, [CII] 158 /im and rotational transitions of CO for temperatures below 3000 
K. These terms are only of importance for z < 1.5. Heating due to photo-electric emission from 
grains is included and contributes for z < 1. At every redshift Equation (9) is used to determine 
the metallicity and expression (8) yields the ambient radiation field. The latter is integrated 
over the hydrogen ionization and H2 photo-dissociation cross sections. The ionization, chemical 
and thermal balance is solved iteratively until convergence is better than 1% in the temperature 
and H2 abundance. The redshift dependence of the baryonic mass fraction is included in the 
calculation of the column density. In deriving the H2 abundances we have included the radial 
N(HI) profiles of proto-galactic disks as presented by Kauffmann (1996). 

5.2. The Moderate Epoch 

The results for a fixed star formation rate S = O.lM yr _1 are presented in Figure 4. Between 
z = 4 and z = 2, the gas is kept at a fairly constant temperature of approximately 500 K. This 
confirms the suggestion by Haiman et al. (1996) that H2 cooling can facilitate the collapse of 
high-mass objects and initiate the formation of structure in the early universe. It is also apparent 
that due to the importance of the stellar radiation field, massive star formation will not occur 
since a cool phase of T < 300 K gas is not present above a redshift of two. In fact, the existence 
of this moderate epoch ensures that the initial cooling time of the disk is not shorter than the 
dynamical time at redshifts larger than 3. If this would not be so, then the star formation rates at 
high redshifts would be much too large. 

These effects are further reflected in the abundance of CO which remains low down to 
a redshift of unity. When the metallicity reaches a threshold of approximately O.O3Z , the 
magnitude of the cooling curve is strongly enhanced. A multi-phase medium is now expected to 
result. Note that to simplify Equations (10) and (11) for the metal enrichment and dust-to-gas 
ratio development we have assumed a Schmidt law with index, p = 1, implying the star formation 
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rate per unit mass, S/M, is a constant. 



5.3. The Phase Transition 

Figure 5 presents phase diagrams for various epochs and star formation rates. For early 
epochs, H2 is the only coolant below 3000 K and only single temperature solutions exist for certain 
pressures. As the metallicity increases the P — n curve attains its characteristic S-shape and lines 
of constant pressure cut it at three different temperatures. That is, a multi-phase medium has 
formed in which the cool and warm components are stable (Field, Goldsmith & Habing 1969) 
and which can facilitate large scale star formation. From Figure 5 one finds that the transition to 
a multi-phase ISM in the bulk of the galactic disks occurs at z « 1.5. At this point it is timely 
to mention the recent observations of galaxies in the Hubble Deep Field presented by Madau et 
al. (1996) and Mobasher et al. (1996) for the star formation history of the universe. Although 
the published results are preliminary, They seem to suggest evidence for a star formation epoch 
around z ~ 1 — 2. 

Assuming that there is 1 supernova per 1OOM of stars formed, the filling factor of the hot 
gas component due to supernovae is (McKee & Ostriker 1977) 

Q = °- 5 ^i 28 ( 1MQ 5 yr _ 1 )( 5 X W M M0 )^)(^F3) L14 ( 2 . 5 x lows*)" 1 ' 7 " ^ 
which can be written using the critical star formation rate rate as 

Q = 0. 5 j^( 5Xl0 " M °)(^r" m 1 '" (^-)-»-« (34) 

* m ttb, g viooy v ioook ; v ; 

For a Schmidt star formation law, and using Equation (1), we find Q oc M e (l + z)^ - 38 where, 
< e < 1/6, and consequently the estimate is fairly robust. 

Therefore, for high star formation rates at the critical value (i.e. more massive objects) the 
filling factor of hot gas can be significant in these protogalactic disks since the pressure is low 
and the supernova bubbles can expand and fill large volumes. The above calculation is indicative 
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only since the star formation is usually clumped and the supernova bubbles generally break out 
of the disk and vent their energy into the halo of the protogalaxy. This process can also result in 
a feedback that inhibits star formation (Norman & Ikeuchi 1989). However high star formation 
rates a multi-phase medium may develop where the hot phase is driven by supernovae energy 
input. For star formation rates of 0.1 — 0.3M Q yr _1 and below, the hot phase is not significant and 
a two-phase mode can occur. 



6. Formation of H2 on Grains 

An additional important ingredient in the regulation of the star formation is the formation of 
H2 on dust grains which renders its abundance independent of the electron fraction. These effects 
have been included in the models above, but an analytic model is presented here to facilitate a 
connection with observations. 

In the presence of dust and at low densities, molecular hydrogen is formed through grain 
surface reactions. In steady state, the local H2 density is given by (Tielens & Hollenbach 1985, 
van Dishoeck & Black 1986) 

»<*>-rftb < 35> 

where nu = n{H) + 2n(i?2) is the total hydrogen density. We have discussed the mean density 
in the disk but it is necessary to make a good estimate of the mean density in the diffuse clouds 
in the cold neutral medium component where the H2 formation is most likely to be initiated. 
We choose to parameterize the dumpiness of the cold neutral clouds by a diffuse cloud covering 
factor f c . If, as we have indicated above, the protostellar disks are not subject to substantial 
massive star formation before molecules are formed then we do not expect the development of a 
supernova-driven hot (~ 10 6 K) component. Consequently, the temperature and density contrast 
between the cloud and intercloud component will be modest. We infer that the covering factor of 
the diffuse cloud component is of order ~ 1. 
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The controlling parameter, Arijj, is the ratio of the formation rate of molecular hydrogen on 
grains to the photo-dissociation rate of H 2 .A is given by 

IuvRmn(3(T)e-^C 1 > 

where lyy denotes the enhancement of the average UV background with respect to the interstellar 
field at z = 0, tuv,C is the optical depth in the UV continuum and tjjv,c = 2. 5£ Ay, Rthin is the 
unattenuated H2 photodissociation rate and (5 is the self-shielding function. The rate constant, 
k s , depends on the nature of the grains. The value of Rthin depends on the precise slope of the 
UV background between 912 and 1100 A. 

For the present discussion, we want to determine for which epoch molecular hydrogen 
self-shields in the cold diffuse clouds created in the phase transition, and therefore adopt the 
constraint 

A > — , (37) 
n H 

yielding a fractional H2 abundance of 0.4. For self-shielding, the constraint (37) reduces to 

IuvRthinT UVtL exp(-T[/ Vt c) < {^jk s T 1/2 ^n H . (38) 
Explicitly including the redshift dependence and the dependence on covering factor we find that 

IuVfiRmnfhuV,L,0 e M-TUV,C,0fcHl + Z)*)(l + Z)^ < Q) ksT^QUHfl. (39) 

Subscripts with zero refer to present time values. 

For S ~ 0.03 — O.lM yr _1 the internally generated radiation field dominates at lower 
redshifts, z ~ 1. In the limiting case where q = 1, the total extinction varies slowly with redshift 
since the lower gas-to-dust ratio of the absorber compensates for its higher column density. In this 
case, the critical redshift for H2 formation, z mo i, is given by 

(1 + Zmol) < {luv fi RiLfc 2 4 v ^ exp (Tuvfififc^ksT^ZonHfi) 2 • (40) 
With a numerical value z mo \ = 1.6. 
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Including variations in the value of S ~ 0.1 — 0.5, this yields a fiducial range for H 2 formation 
of, 1.2 ^ z <; 2.0. These estimates have the advantage that they are weakly dependent on the input 
parameters. The strongest dependence is on covering factor and star formation rate. Obviously, 
for small covering factors molecules may form at higher redshift but they will be less observable. 
Note that the epoch of significant dust shielding and cooling are similar since the dust production 
and metallicity production are closely related. 

The increasing dust columns can absorb a significant fraction of the internal stellar radiation 
field and this will influence the the time scale on which subsequent star formation proceeds. Using 
the formalism of McKee (1983) it is possible to estimate the effect of the dust column on the star 
formation rate. There is a critical extinction of approximate 4 mag, necessary to absorb most of 
the radiation of newly formed stars. The time scale to convert all available mass into stars for a 
density of roughly 1000 cm~ 3 is given by 

3.2 x 10 7 a^)°- 5+1 - 25/Av [(1 + ^)°- 5 + (^)0-5]l-2.5/>lv e 16Mv yr (41) 

and the dust extinction Ay = 4£(1 + z) 2 . For z > 2 this time scale is much larger than the local 
Hubble time, but for epochs later than 1.7 star formation can proceed efficiently in the bulk of the 
galaxies. 



6.1. Statistical Approach 

The covering factor discussed above is ultimately determined by the process of galaxy 
formation. Using the statistical approach discussed in Section 2, for any given redshift, the 
condition (37) for H2 formation on grains can be written as 

r > r ,crit(z). (42) 
The fraction of galaxies which contain molecules, f mo \(z), is then given by 



Jmol(Z) = 755 Z^Tl 1 (46) 

Je(o,z) e ax 
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where £(r ,z) = ( 2^a{r ) ) ' definition of the statistical epoch of molecule formation is taken 
to be /moi(^) ~ 0.5. Inversion of this condition yields < z mo \ >. For the parameter values adopted 
above and putting the biasing parameter 6=1, the relations ro jC rit = 0.17(1 + z) 2+2a Mpc and 
< z moi >= 1-5 approximately hold. Note that ro iCr it does not dependent on biasing, but the value 
of < z mo \ > does. 

Finally, the probability, V, that a given line of sight to a QSO will show molecular absorption 
by an intervening damped Lyman Alpha cloud in a (proto-) galaxy is given by (Peebles 1993) 

V » J anocH^n' 112 ^ + z) l ' 2 dz, (44) 

where the integral extends over redshifts between and < z mo \ > and a = -irr 2 .. The characteristic 
radius (bright part) of the galaxy is chosen equal to r g = 10/i _1 kpc. With the parameters given 
above this yields V « 8 x 1CT 2 . 

Once the conditions are satisfied for molecular hydrogen formation the time scale is of order 
~ l/(nH^)Gyr ~ 10 8 (1 + z)~( p+1 ^yr. The ion-chemistry proceeds much faster on a time scale of 
approximately 3 x 10 3 (1 + z)~ e yr. Thus, once the hydrogen molecules can be formed on grains, 
the full range of diffuse cloud species follows immediately. 



7. Conclusions, Discussion and Implications 

Initially motivated by QSO absorption-line observations indicating an absence of molecular 
hydrogen at high redshift, we have examined the formation of molecular hydrogen, after the epoch 
of re-ionization, in the context of canonical galaxy formation theory due to hierarchical clustering. 
The issue of molecular hydrogen at high redshift has been discussed in recent interesting papers 
by Haiman et al. (1996) and Tegmark et al. (1996). Whereas they concentrated more on the 
redshifts above z ^ 10, we are focusing on the physics of the phases of the protogalactic gas at 
reshifts z ~ 1 — 10. We have been able to use recently computed accurate models for the cosmic 
background radiation field (Haardt & Madau 1996) and for the radiation field from a stellar 
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population with a Salpeter IMF and with a 1M© lower-mass cut-off in a continuous star forming 
mode (Leitherer & Heckman 1995). We have also used a simple model for a protogalactic disk 
and followed conventional CDM cosmology extended by more recent numerical calculations when 
considering the spectrum of disk masses and column densities. We have mainly concentrated on 
structures associated with la and 2a density fluctuations which are expected to form the bulk of 
the stars. 

As with the recent work by Haiman et al. (1996) and Tegmark et al. (1996) we find there 
is an initial epoch of H2 production in the gas phase produced through the H~ channel route 
where the abundance of molecular hydrogen is approximately 1% and given simply in terms of 
the ionization parameter by equations (20) and (22). Predicting the details of the state of the 
interstellar medium of protogalaxies is a complex task. We have normalized to what is known 
from work on the neutral phase of our Galaxy and assume that the fundamental units in the gas 
phase of protogalaxies during this epoch are akin to the diffuse clouds found in our own galaxy. 
We have shown that until the metallicity of the gas achieves Z ~ 0.03 — 0.1 at a redshift of 
1-2, cold giant molecular clouds are not formed due to inefficient cooling. We have found that 
star formation in the protogalaxies can become self-regulated due to heating of the gas by the 
internal stellar radiation field. We have given a simple analytic model for the feedback process in 
Section 4.1. It is possible to define a maximum star formation rate during this epoch. Plausible 
estimates give a rate of <^ lM Q yr _1 and therefore we have considered an appropriate fiducial star 
formation rate to be ~ 0.1 — 02M & yr l . The production of metals and dust proceeds slowly in 
this self-regulated mode. This slow star formation phase was shown to terminate once the metal 
abundance increased to a level of approximately Z ~ 0.03 — O.IZq. From an analysis of the phase 
diagrams in Figure 5, we found that species such as C, O and CO become sufficiently abundant 
and can cool the gas below 300K to ~ 10 — 30K. At this point a phase transition can occur 
in the protogalactic gas. For the low fiducial star formation rates discussed above, we find this 
to be a transition to a two-phase medium as described by Field, Goldsmith & Habing (1969). 
Dense molecular clouds can form and the star formation is no longer self-regulated in the manner 
described above since the UV radiation does not penetrate the dense cores of the clouds. We 
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expect that rapid, massive star formation ensues and the abundance of metals and dust increase 
concomitantly. The dust abundance also becomes sufficiently high to allow molecular hydrogen 
formation on grain surfaces. With the increased star formation rates, the ISM will change to one 
dominated by supernovae energy input (McKee & Ostriker 1977) with significant exchange of 
mass, energy and metallicity from the disk to the halo (Norman & Ikeuchi 1989). In a subsequent 
paper we will investigate in more detail the effects of this phase transition on the evolution of 
dwarf galaxies and the importance of metal loss driven by supernova explosions. 

Our analysis may be of relevance to the G-dwarf problem: the Simple model (Pagel 1989, and 
references therein) for chemical enrichment overproduces the number of metal-poor stars (Cowley 
1995; Worthey, Dorman & Jones 1996 and references therein). Of the many solutions proposed 
for the G-dwarf problem, the simplest appears to be that by the time a few percent of the gas 
mass of a galaxy is assembled into stars, the remaining gas reservoir is already enriched in metals 
and has not yet experienced any star formation. The moderate phase due to feedback identified 
in this work, is likely to cause star formation in an inhomogeneous way. The very nature of the 
feedback mechanism dictates that star formation in one location strongly suppresses additional 
star formation in its vicinity. The subsequent phase transition then causes rapid star formation 
throughout the mostly unprocessed ISM which now has a metallicity close to 10% of solar. 

In our analysis we have assumed a fixed value for the collapse factor A -1 . In reality, the 
distribution of spin parameters may be quite wide (Warren et al. 1992; Dubinsky h, Carlberg 
1991). If we view the disk as being formed from a spherical object sustaining a low H2 abundance, 
driven by the extragalactic UV background, then an increase in the collapse factor will increase 
the total and H2 column densities quadratically and the local density like A 3 . Consequently, more 
of the stellar radiation can be absorbed locally, preserving the H2 abundance and suppressing 
the global heating, and a higher star formation rate can be sustained. We estimate, although 
tentatively, that for objects with A ~ 0.02, a factor of 3.5 below our canonical value, the redshift 
at which the phase transition occurs, can be as high as three. More detailed knowledge of the A 
distribution is necessary to assess how common such high redshift objects are. 
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In summary, from our elementary cosmological model we conclude that this new mode of star 
formation, where objects now akin to giant molecular clouds in our Galaxy become the sites of star 
formation, occurs at a redshift of approximately 1.5 with a value higher by a factor of 2 if more 
massive initial perturbations or larger collapse factors are considered. The phase transition in the 
interstellar medium of protogalactic disks as analyzed in this paper is now a plausible physical 
reason that the formation of disks of galaxies occurs at a redshift of order unity with a significant 
increase in star formation after the metallicity has achieved a value of order Z ~ 0.03 — 0.1Z Q . 
These findings are consistent with the recent studies of the Hubble Deep Field (cf. Madau et 
al. 1996; Mobasher et al. 1996). The combination of feedback and a phase transition can provide 
a natural solution to the G-dwarf problem. 
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populations and the cosmic background radiation field. We are also grateful to Piero Rosati for 
his assistance in the presentation of the numerical results. We thank the anonymous referee for 
his detailed and valuable comments. MS also acknowledges with gratitude the support of NASA 
grant NAGW-3147 from the Long Term Space Astrophysics Research Program. 
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Fig. 1. — Figure 1 shows the abundance of halos as a function of mass and redshift for standard 
CDM with biasing parameter 6=1. The lower M — z plane shows, as functions of redshift, the 
characteristic turnover mass, (solid curve), in CDM and the limiting mass of the dark halo, 
M coo i (dashed curve), in which gas can cool. 

Fig. 2. — Figure 2 shows the combined stellar and background radiation fields at redshift z = 1, 
2.5, and 4 for star formation rates of 0.01, 0.1 and lM yr -1 . For rates larger than 0.03 the stellar 
component dominates and feedback is important. 

Fig. 3. — Figure 3 shows the derived gas pressure, radiation pressure and ionization parameter as 
functions of redshift for S = O.lM yr -1 . The radiation pressure is typically much less than the 
gas pressure due to the low star formation rate at which the feedback mechanism is most efficient. 

Fig. 4. — Figure 4 shows the results of the calculations for the gas temperature and the abundance 
of molecular hydrogen, CO, dust (£ g d), and metals (Z m ) as functions of redshift for a star formation 
rate S = O.lM yr _1 . 

Fig. 5. — Figure 5 shows pressure, P, and gas density, uh, phase diagrams for redshifts 1, 2 and 3, 
respectively and star formation rates of 0.01, 0.1 and lM yr _1 . The curve in the upper left panel 
denotes constant pressure. The symbol Z indicates the metallicity with respect to galactic at that 
epoch for a particular star formation rate. The labels F, G and H denote the thermal equilibria of 
which G is unstable. 
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